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The global solar magnetic field — identification of travelling, long-lived ripples 

R. K. Ulrich and Tham Tran 
Department of Physics and Astronomy, University of California, Los Angeles, CA 90095 

ABSTRACT 

We have examined the global structure of the solar magnetic field using data from 
the Fel spectral line at A5250.2A obtained at the 150-foot tower telescope at the Mt. 
Wilson Observatory (MWO). For each point on the solar surface, we find the value of the 
magnetic field in the meridional plane, i?m, by averaging over all available observations 
using a cosine weighting method. We have revised our cosine weighting method by 
now taking into account more fully the highest latitude geometry. We use the annual 
variation in the latitude of the disk center, bo, to deduce the tilt angle of the field 
relative to the local vertical so that we can find the radial component of the field, B^ 
from Bm- We find this tilt angle to be small except for a near-polar zone where a 
tilt-angle model can reduce the annual variation. The reduced annual variation in the 
deduced B^ allows us to study dB^/dt and associated deviations in B^ from a smoothed 
Br with a smoothing width of 2.5 years. These functions make evident the presence 
of small amplitude (3 to 5 gauss) but spatially coherent ripples with a semi-regular 
periodicity of 1 to 3 years. At any given time the half-wavelength (peak-to-trough) is 
between 15 and 30 degrees of latitude. These patterns are ubiquitous and in many cases 
drift from near the equator to the poles over a time period of roughly 2 years. The drift 
rate pattern is not compatible with simple advection. 
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1. Introduction 



Magnetic and Doppler measurements of the full solar disk have been obtained and recorded 
digitally at the 150-foot solar tower telescope at the Mt. Wilson Observatory since 1967 (hereinafter 
the MWO synoptic program, note that earlier observations were recorded by means of various analog 
systems but these have not been incorporate d into t he digital record). Most a nalyses of the s e data 
such as those given by ISnodgrass fc Ulrichl (|l99d ): lUlrich fc BovdenI (|2005l ): lUlrich et al.l (|l988l . 
2003) have relied on a quick-look style database we call the Interactive Data Reduction (IDR) 
which includes r nagneti c and Doppler arrays having spatial dimensions of 34 by 34. The recent 
analysis given by lUlrichl ()20ld ) was based on a new approach which starts from daily average arrays 



having spatial dimensions of 256 by 256 wherein ea ch pixel is co rrected for differential rotation prior 
to its inclusion in the average. We also include the lTranI (|2009l ) extrapolation for pixels in the polar 
regions that are not imaged due to the inclination of the sun's axis of rotation to the line-of-sight: 
60 • In addition we have used a superposed epoch analysis to find a model for the geometry of the 
solar surface magnetic field which minimizes the annual variations due to the changing 69 • Because 
of these improvements, we are able to study small amplitude variations in the radial component of 
magnetic field Br- We report here on a new phenomenon that emerges from this cleaned- up data 
set: a ubiquitous set of large-scale wave-like features ha ying cycle peri o ds ra nging between 0.8 and 
2 years. Similar features have been found recently by IVecchio et al.l (|2012l ) in the magnetogram 
data from the US National Solar Observatory. These features are compared briefly with previously 
reported periodicities in solar phenomena, called solar Quasi-Bienniel Oscillations (solar QBO's), 
but are not studied in enough detail to reveal any physical relationships. 

We start with a summary of the magnetic field data and the methods of correction and aver- 
aging we have used. The following section presents the wavelike features and shows that the field 
patterns are not simply advected relics of active region. The final section offers a few thoughts 
about what these features might represent. 



2. The Magnetic Field Data 
2.1. The initial data 



The Babcock magnetograph at the 150-foot tower telescope obtains full-disk magnetograms by 
scanning the solar image over the entrance slit and after passage through the spectrograph samples 
24 spectral channels 400 times per second in two states of circular polarization. Averages of the 48 
numbers are obtained while the solar image is moving over the spectrograph entrance slit and stored 
in a series of records each of which includes the 48 numbers along with the x, y coordinates of the 
entrance slit on the solar image as well as the average time for the observed number. The observed 
data are subsequently assembled into a rectangular grid of pix els of 256 by 256 a t a pa rticular time 
using the differential rotation correction method described bv lUlrich &: BovdenI (J2006l ). 



The raw data numbers obtained from the MWO system are initially stored as 48 polarized 
intensities along with the red and blue stage Doppler servo positions as functions of x, y, t where x 
and y are the horizontal and vertical positions of the entrance slit on the stage and t is the average 
time for the pixel. We concentrate here on the case of the line at A5250.2A and use the two spectral 
samples on opposite sides of that line so that each pixel yields four quantities - two intensities for 
each spectral sample. The four intensities are corrected for blue and red stage servo error and 
converted according to center-to-limb dependent calibratio n curves into a magn etic signal and a 



Doppl er signal. These steps are described in greater detail in [Howard et al.l (jl983l ) and 
(I2OO2I ). Our applications of the d ata to modeli ng of the Total Solar Irradiance (TSI) ( 



Ulrich et al. 



Ulrich et al. 



2OIOI ) and meridional circulation (jlJlrichj |2010J) have demonstrated the effectiveness of using daily 



averages and limiting the days used to those days having three or more observations and we have 
adopted that approach in the present analysis. Prior to 1986 there were few days that satisfied this 
condition and all good observations were included in the analysis. 

The observations are made using one of two entrance aperture, slow-grams with a 12 arc-sec 
squared aperture and fast-grams with a 20 arc-sec squared aperture. Comparison of the magne- 
tograms shows that the fast-grams have a systematically smaller field than the slow-grams. Our 
program has a mix of slow- and fast-grams and we find that the systematic effect can be removed 
by dividing the field from the daily average (-Bios)ave by a factor Cfs to get an equivalent field that 
would have been measured as a slow-gram (-Bios)siow We find that Qg ~ 1 — 0.0127 A'fast where 
-^fast is the number of fast-grams on that day. Typically A^fast = -^totai — 2 where A^totai is the 
number of observations of both types durin g the day. In a dditio n, we have multiplied the observed 
magnetic fields by the correction factor by lUlrich et al.l (120091 ) to bring the field strength to the 
value appropriate for the base of the chromosphere. As discussed in the above reference this factor 
corrects for effects due to the fluxtube nature of the photospheric magnetic field. 

The array size we use for this analysis is 256x256 yielding a pixel size of about 8 arc-sec. This 
is smaller than the actual pixel size but due to the dithering effect of adding together multiple 
observations, this pixel size is close to an effective resolution element. We have assembled various 
results from the analysis into latitude dependent arrays by averaging into bins 2° wide with the 
bins being centered on odd values of latitude (this avoids dividing by the sine of the latitude at 
the equator and by the cosine of the latitude at the pole). We will generally use degrees instead of 
radians to discuss the dependency of Br on latitude. 



2.2. Calculation of the radial field from the line-of-sight field 

The de duction of the magnet ic field geometry h a -s bee n the topic of a number s tudies such 
as those by ISvakaard etaP hdjj ). IWang fc SheelevI (|l992l ) and IPetrie fc Patrikeeval (I2OO9I '). We 
describe here a new method that depends on the long duration of the di gital record from the 
MWO synoptic program. Following an approach similar to that described by lShrauner &: Scherrer 
(JI994I ) we use solar rotation to resolve the line-of-sight magnetic field B^ into components in the 
meridional plane B^ (the plane through the observed point and the solar axis of rotation) and 
in the zonal plane B^, (the pla,ne par allel to the solar equator). Our method is described most 
completely bv lUlrich &: BoydenI (120061 ) and results from our method have been presented at various 
times starting with lUlrichI ( 19931 ). The essential ingredient in our approach and the Shrauner- 
Scherrer method is to weight the observed line-of-sight field Bi by the sine and cosine of the central 
meridian angle 6L and sum over all available observations. The meridional and zonal components 
then can be recovered from a matrix inversion. We have found that the inclusion of sums of cross 
terms involving the products cos{SL) sm{SL) magnifies errors in cases where the data is sparse. For 
uniformly distributed data, these terms should vanish but in our actual data sets, they do not. We 
find the final result is improved by setting such cross terms equal to zero. 



In this paper we are interested in the magnetic field at high latitudes and a term neglected 
previously needs to be included. Equations relating magnetic field components in a rectangular 
gri d for an observe d image to the components in a heliographic coordinate system have been given 
by iHagvardl (|l987l ) . Because our observations are made after we have determined the east- west 
direction at the focal plane and oriented the scan lines so they are perpendicular to the sun's axis 
of rotation, our images always have the p angle effectively set to zero. Also we only measure the 
line-of-sight field component Bi. The Hagyard equations can then be written: 



Bf 



[{Br COS b + Bg sin b) cos 6L — B^ sin 6L] cos 60 + {Br sin b — Bq cos b) sin 60 (1) 



where b is the latitude, bo is the latitude at disk center, Br, Bg and B^ are respectively the radial 
component of the magnetic field, the component of the field in the colatitude direction and the 
component of the field in the (j) direction for a right-handed spherical coordinate system. We define 
a tilt angle (" = arctan{Bg/Br) giving the angle of the field line with respect to the local vertical 
measured positive from the north polar axis toward the south polar axis. We describe below how a 
model for ^ as a function of b can be found by minimizing the variance of Br- Using ^ equation([T]) 
becomes: 



Bi cos C = [Br cos (6 — C) cos 6L — B^ cos C, sin 5L\ cos 60 -|- Br s\n{b — C) sin 69 



(2) 



The next step is to multiply by cos5L and sum over all available observations. In accordance with 
our modification of the Shrauner-Scherrer method, we set ^cos5Lsm.6L = and ^sva.6L = 
and assume that Br is constant over the set of observations. This gives: 

cos CJ2^e cos 6L 



Br 



cos(5 — C) cos 60 l^(cos 6L)'^ + sin(6 — () sin 69 Yl cos 5L 



(3) 



2.3. Filling in unseen portions of the polar regions 

A new step in our reduction is the filling of th e unseen parts of the polar regions with a 
method based on the ideas developed by iTranI (J2009l ) . The problem of estimating the magnetic 
field for those portions of the disk that are beyond the limb has been addressed in a variety of ways 
and having a good approach is important for the polar regions whose dipole-like field influences a 
major portion of interplanetary s pace. A var i ety of methods have been developed and tested with 
an excellent recent discussion by ISun et al.l (|201ll ) which includes a summary of seven methods 
and introduces a new extension based on a combination of the best features of several of those 
methods. The method we present here differs from all of those by transformin g the i mage to a 
60 = configuration prior to the step of deriving the fill values. The study by iTranI (J2009l ) has 
shown that this approach avoids generating field structures aligned along lines of constant longitude 
which are seen in some polar fills based on polynomial interpolation. The Tran method as we have 
implemented it is stated below; and we also illustrate the performance of the method using samples 
of individual maps as viewed from the pole. A detailed comparison between this method and other 
methods is, nonetheless, beyond the scope of this paper. 
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Fig. 1. — This figure shows a sequence of projections of maps of the radial component of the sun's 
magnetic field. Each circle shows By as a function of position for a Carrington time value as indicated 
below the circle. The Carrington rotation number and corresponding Carrington longitude apply to the 
line extending vertically down from the pole point to the lower edge of the circle. For other angles the 
heliographic longitude can be computed from angle difference from this vertical line. This longitude is not 
the Carrington longitude for any other angle on the circle. All points on the fig ure have been corrected t o 
correspond to the labeled Carrington time according to the methods given bv lUlrich &: BovdenI (J2006l ). 
The latitude as projected is given above the lower right circle and is shown on this figure by the dashed 
circles. Each circle is confined to latitudes pole- ward of 60°. The circular solid lines show the boundary of 
the zone where the field value must be determined from our polar fill technique. Each circle is offset from 
the previous by an interval that corresponds to the rotation period at the pole. The curve on the lower 
right circle gives the pattern that would be produced on that circle if a vertical straight line descending 
from the pole point had been drawn on the preceding circle. This curve can be used to estimate which 
features might have been preserved from one polar rotation to the next. 



The Trail method is particularly well suited for application to the MW O data because 1) the 



error of measurement from the MWO system is substantially less than 1 G (jUlrich et al.ll2009l ). 2) 
the scanning aperture method of measurement prevents scattered light within the spectral analyzer 
from influencing the field strength near the limb in contrast to filter-based system where scattering 
from bright parts of the image can reach the limb, and 3) the signal is sampled at 400 Hz so that the 
seeing fluctuations are cancelled in the algorithm that provides the magnetic signal. An empirical 
estimate of the noise is difficult and our best resource is a set of images taken periodically with 
the circular polarization analyzer turned off. The distribution of apparent magnetic signal from 
these pseudo-magnetograms gives a gaussian width of 0.5 G. Each daily average is composed of an 
average of 1 1 individual magne tograms while near the limb, the fluxtube physics correction factor 



discussed in lUlrich et all (|2009l ) multiplies the raw field strength by about 1.8 and the conversion 
from line-of-sight to radial multiplies the field by another factor of 3. The combination of these 
factors gives an error of about 1 G in Bj- for each pixel on the image near the limb. 

To calculate fill values we take the observed values of Br derived from the geometry described 
above and rotate the image so that the equator is at the center of the disk. The observed pixel values 
then occupy positions on this image different from where they were observed. In particular pixels 
near the pole do not have observed values and these are the pixels for which we seek appropriate 
estimated values. For all such pixels, we find on this rotated plane image the six nearest pixels 
with observed values. We then calculate the projected distance on this rotated image between 
each observed pixel and the pixel for which we wish to find an estimate. We take for this estimate 
the weighted average of the observed pixels with the weighting function being the inverse distance 
between the target pixel and each observed pixel on the projected and rotated image. This method 
allows longitudinal structures to influence the nearby polar pixels in a natural and bounded manner. 
In the 6o = geometry, the nearest pixels as defined this way include a range in longitude which 
increases as the pixel is further from the edge of the observed image so that there is less longitudinal 
concentration of the nearby structures. 

In order to illustrate the performance of the pole filling approach, we have selected times near 
polar field reversal to form eight snapshot maps of the radial component of the solar field. The 
coordinates on these maps are heliographic longitude and latitude. Each pixel is corrected for 
differential rotation and averaged according to equation ([3]). The map has been reprojected so that 
it appears as if being viewed from above the sun's north pole. Rotation is in a counter-clockwise 
direction and the time of all points is corrected to that on the central meridian (a line vertically 
downward from the circle center). Each circle shown only includes l atitudes northwa rd of 60°. 



The format of the figures have been chosen to agree with that used by ISun et al.l (|201ll ) although 



the times chosen do not corresp ond to those used by t hose authors. The time of observation is 



given in Carrington time units (jUlrich &: BovdenI |2006| ) and is listed above or below each circle 



The corresponding times in fractional years are 2000.849, 2000.952, 2001.056, 2001.160, 2001.264, 
2001.367, 2001.471 and 2001.574. 

We see from the sequence in figure ([T|) that the structures in the polar region during this stage 



of the cycle are quite variable from one rotation to the next. Our method of calculating B^ as given 
in equation ([3|) can artificially enhance features in some longitude ranges as a result of irregularly 
space observations. Some of the strong patterns in the upper left and upper right circles may be 
due to this effect. However, the extension into the unseen areas in the upper right circle is not 
enhanced beyond the features in the seen areas and no fill has been applied to the upper left figure. 
In general the fill pattern is a logical extension of the seen patterns with no evident enhancement 
of noise patterns. 



2.4. Determination of the tilt angle ( 

For an analysis in which (^ is assumed zero (the field lines are exactly vertical everywhere), we 
find there to be a strong annual variation in the derived B,., especially in the polar regions where 
there are long periods during which one magnetic polarity or the other dominates. The influence 
of bo comes from the denominator of equation ([3]) which for 5L near zero is just cos(6 — C — bo). 
As b — (^ — bo approaches ±90°, the line-of-sight magnetic field is multiplied by a factor that can 
be quite large. We isolate this effect by considering an annual superposed epoch analysis whereby 
all the derived values of \Br\ are treated as functions of the time of the year and averaged over the 
full data set. To illustrate the need for inclusion of the tilt angle C we show in figure ^ a sample 
superposed epoch result for latitudes b = 85°N/S. The top panel shows the result for vertical fields 
with C = everywhere. At the latitude chosen, cos{b — bo) ~ sin(6o)- For comparison we show bo 
in this top panel with the scale adjusted in amplitude and offset to coincide with the superposed 
epoch result for -Braverage ^t the south latitude. 

We can explain the nature of the pattern in the top panel of figure ([2]) by concentrating on 
the plot for the southern hemisphere where there is greater regularity. The sense of the variation 
is such that when bo is the most negative and 85° S is most visible, the deduced field is relatively 
weak. When the latitude is least visible and is filled in by regions slightly to the north, the result is 
stronger. On the adjacent latitude of 77°S where filling by extrapolation is not needed, the result 
is a relatively constant 6.5G without the 2.5G annual variation seen at 85° S. Interpreted literally 
this behavior implies that there is a decrease in the field strength at the pole and that the field 
strength gradient depends on the time of the year. We believe that the correct interpretation is 
that the field is not radially oriented near the pole and that the correct from Bg to Br requires 
a larger multiplicative factor - hence the field is tilted slightly away from the line-of-sight and is 
more nearly perpendicular to the line-of-sight. We impose an assumption that the field lines do 
not converge toward each other; but, if some model were to suggest that this assumption is invalid, 
it would not impact the effectiveness of our analysis. This assumption imposes the condition that 
C, > b — 90° for 6 > and C, <b + 90° for 6 < 0. As a convenience in developing our model, we have 
used the quantity ^ = 90° — |6| to represent the proximity to the poles. The model for C, we have 
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Fig. 2. — This figure shows a superposed epoch analysis of Br as derived from equation ([3]) and 
averaged over the full data set. The superposed epoch average is obtained by adopting the inde- 
pendent variable to be the time modulo 1 year. The portion of the figure between 1.0 and 1.2 year 
is a copy of the function between 0.0 and 0.2 and is included to facilitate examination of the trends. 
The top panel has the tilt angle C taken to be zero so that the assumed field orientation is exactly 
vertical. A selected pair of latitudes at 85°N and 85°S are shown along with a plot of the apparent 
polar tilt angle 6o as the solid line. The scale of 6o and location of the line has been adjusted to 
aid in the comparison to the case of 6 = 85° S. The lower panel shows the two latitude results after 
a model for C as described in the text is applied. 



adopted is: 

C = sign(5)max[Coexp(-(eM)'),l^|-90°] (4) 

where sign = 1 for 5 > and sign = — 1 for 6 < and where max[a, b] is the maximum of a and 
b. We have adopted Co = —6° and 5^ = 12.5°. Figure ([3]) ihustrates the geometry of the field lines 
near the sun's north pole. The improvement in the variance of \Br\ is illustrated in figure dH). 

2.5. Relationship with the fiuxtube correction factor 

As indicated above, the analysis presented here has been applied to the line-of-sight magnetic 
fields that have been cor rected for the effec ts of photospheric fiuxtube physics according to the 



recommendation given bv lUlrich et al.l (J2009l ). In keeping with the traditional approach with MWO 
data, we initially worked with the uncorrected fields. With that approach, we found ( to be near 
zero close to the poles and to have a maximum of about 13° for latitudes near 65°N/S. We were 
surprised to find that after the application of the fiuxtube correction factor, the field became radial 
almost everywhere with the exception being the very high latitude regions presented above. We 
checked to determine if the effect could be due to an error in the treatment of the near-limb pixels 
around the whole circumference of the solar image. This last assumption failed to remove the near- 
pole effects because the feature we correct is actually a contrast between properties of the solar 
surface at latitudes greater than 80° compared to those between 65° and 80° and those latitudes are 
altered equally in the case of a full-circumference correction. The character of the center-to-limb 
dependence of the field strength is discussed further in the subsection below. 



2.6. Influence of the center-to-limb angle 

Although the infiuence of the bo angle on the near-polar field is most strongly seen in the polar 
fill zone, in fact the dependence of the fully seen pixels on both the latitude b and the center-to-limb 
angle p illustrates the need for the Co = ~6° model we have adopted. Due to the long temporal 
record available to us, we can study the systematic dependence of the field on both these variables 
by collecting all observed points into bins of b and p. We have taken the bin width to be 2 degrees 
in extent centered on odd values in degrees. 

The simplest approach is to average the absolute values of the Br but this approach produces 
an artificial increase in | -Br I average due to the noisiness of the points nearest to the limb where the 
noise can be comparable to the value of Bj.. We can overcome this effect by using the fact that in the 
near-polar zone of 70° to 90° and in the high-latitude zone of 50° to 70° , the field is mostly positive 
or negative between the times of polar field reversal. Thus we can average the signed value of -Br- 
and obtain nearly the same result as averaging the absolute value as long as we reverse the sign 
whenever the polar region has a negative predominant polarity. Then when the noise is comparable 
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Fig. 3. — This figure shows the geometry of the adopted field line tilt as defined by the tilt angle 
Q. Each of the arrows is in the direction of the field lines but for purposes of illustration are shown 
with a length of 5% of the solar radius. For reference, the longer arrows show the direction of the 
local vertical. Each of the shorter arrows is at the center of one of the bins used to describe the 
magnetic field pattern. The C character with the curved arrow points to a place where the tilt angle 
has the largest value. The x and y coordinates are in units of the apparent solar radius and are 
oriented in the same way as the image is observed. The expanded inset shows the definition of C,. 
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Fig. 4. — This figure shows the variance of the superposed epoch averages of \Br\ (see figure[2]) as a 
function of latitude in the polar regions for the case where the field is assumed vertical everywhere 
and for the case where the field obeys the model given by equation Q. 
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Fig. 5. — This figure shows the latitude dependence of the latitude restricted center-to-limb function 
as defined in the text. The points averaged together are all observed and not a result of the fill 
procedure. The error bars are the errors of the mean derived from the variance of the data points 
included for each point. The two lines show averages from deduced values of B^ based on two 
models of the field geometry: the first with a strictly radial assumption (the dashed line) and the 
second with the poleward tilt implied by the model discussed in the text (the solid line). 



12 



to the value, we will not get an increase in the average because the positive and negative offsets 
due to noise will average to zero. The dates we have adopted are: prior to 1981.0 the sign in the 
south is reversed, between 1980.4 and 1990.2 the sign in the north is reversed, between 1991.5 and 
2000.7 the sign in the south is reversed and after 2000.7 the sign in the north is reversed. These 
times are appropriate for 70°. Because the fields are small near these reversal times, the shifts in 
time of reversal for other nearby latitudes do not impact the averages. In addition, after average 
is computed for both north and south latitudes, we average together results for equal values of \b\ 
and call the result the Latitude Restricted Center-to-Limb Function (LRCtLF): {Br{\b\, p)) . The 
restriction in latitude also restricts the range of possible values for p since the closest any latitude 
gets to the disk center at any time of year has p > \b\ —bo. 

Because the LRCtLF isolates observations having the same values of p, the resulting latitude 
dependence can only come from the combination of magnetic field strength and magnetic field 
geometry. While we cannot separate these two factors, we can compare the implied functions based 
on the C model. In particular we examine the cases Co = 0° and ^o = —6°. Figure ([5]) shows that 
the strictly radial field assumption results in a noticeable reduction in field strength at the highest 
latitudes whereas the model with the poleward tilt shows a pattern of field strength which increases 
smoothly toward the pole. 

For the highest latitudes, the LRCtLF is difficult to use to study the p dependence because p 
can only have values between 90° and \b\ — 7.25° so the functional pattern of the dependence of Br 
on p cannot be studied. In addition we have too few points for p > 82° to provide a statistically 
valid result. The variance test of the superposed epoch values of \Br\ includes the result of the fill 
technique and brings in information from lower latitudes. Over the restricted range of p available 
for the observed points as opposed to the result including the fill points, we find that from lowest 
p to highest p, the LRCtLF decreases by 11% with (q = —6° and 23% with (^o = 0°. This smaller 
change in B^ with p is responsible for the lower variance in the superposed epoch function. We 
have tried modifications to the p dependence of the field strength similar to those applied above 
in section 12.51 and found that changes causing (q to become near zero also result in the need for 
positive C in latitudes below the near-polar zone. We cannot rule out such models but feel the 
present model is simplest in that it requires only one adjustment - the near-polar (" model instead 
of two - an adjustment to the fluxtube physics correction function and a lower latitude model for 
C- We note, however, that the profile for the Fel line at A5250.2A has not been measured with 
our system with p above 70° so that some adjustment to our result could be necessary based on 
improved observations of this line profile. 



3. The Moving Magnetic Patterns 

Long-term trends in the solar magnetic field are usually shown in 2-d plots giving the field 
strength as a function of time and latitude. The MWO project has published a number of these as 
derived from the IDR records whose basic spatial resolution starts with a grid 34x34 in sin(central- 



13 



Carrington Rotation 
1650 1700 1750 1800 1850 1900 1950 2000 2050 2100 




1980 



1990 



Year 



2000 



2010 



Fig. 6. — This figure shows the strength of the radial component of the magnetic field as a function 
of latitude and time. Time is given in years along the lower axis and in Carrington rotation number 
along the top axis. The color-bar to the right shows the correspondence between the colors and 
the strength of the radial component of the magnetic field with the units of the field (flux density) 
being in Gauss. Three gaps in the data due to weather conditions and equipment failures have 
been filled by means of a cubic spline interpolation at locations indicated by the grey rectangles 
near 1990.2, 2003.5 and 2010. 
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meridian angle) and sin (latitude). An example is found in lUlrich fc BovdenI (J2005l ). As described 
above, our current reduction differs from those previously as a result of starting with a spatial 
resolution of 256x256 as derived from daily averages of all images available. The prior approach 
yields arrays in latitude that are chosen so that the spacing is roughly uniform in sin(6). That 
produces a coarse result at high latitudes where the bins at indices 1 and 34 include areas between 
the poles and 76°N/S. In addition, the prior plots show the raw, line-of-sight magnetic field instead 
of the radial field after the correction for flux-tube effects. Two consequences of this change are 
1) the field in the polar regions is now significantly stronger than in earlier plots and 2) the 
annual effects seen before are now largely absent. Furthermore, our adoption of the Tran polar 
extrapolation method provides an improved level of detail in the time dependence of the polar 
fields. Our new map of the radial component of the magnetic field flux density Br as a function of 
latitude and time is given in figure ([6]). 

The new map is largely free of the annual variations produced from the variation of Bq, however 
there are a few features worth comment. First, there are sequences like those near the north pole 
between 1992 and 1995 where a four successive local minima are separated by one year. However, 
this sequence and others like it do not persist and do not influence the overall superposed epoch 
values. Second, 1 to 2 years prior to each cycle's final polar field reversal there is a preceeding 
mini-reversal that does not remain. These are consistent in time relative to the reversal even for 
cases where the north and south reversals are themselves out of phase. The color mapping we use 
makes these features more evident because we make the color gradient greater near zero field in 
order to better display weaker fields in the mid-latitudes. In the polar regions when the field is 
reversing, this color mapping combined with the magnetic field ripples discussed below produces 
an effect which appears to be a mini-reversal preceeding the main reversal by one to two years. 

An important property of the Br map is the regular occurrence of magnetic plumes that start 
near the active regions or near the equator and migrate poleward. These are easily recognized in 
maps of this sort and play a critical role in the reversal of the polar magnetic field near sunspot 
maximum. What has been less evident is their ubiquitous presence at all phases of the solar cycle. 
In areas of weaker field, the regular, shorter-term ripples would be lost without the consistent 
quality of figure ([6]). These ripples can be seen more clearly by examining a subset of line plots 
of Br vs. t as are shown in figure ([7j). Although the 2° bin size oversamples the MWO images, 
the consistent appearance and time drift associated with the features indicates that they are not 
simply random noise. 

The ripples can be displayed effectively by examining dBr/dt since this operation enhances the 
short-term changes relative to the 11-year solar cycle pattern. Initially we were searching for times 
of rapid field change in order to find patterns in the Doppler velocities. We believed that the small 
variations in the field were in fact noise so that we reduced these by smoothing. Subsequently we 
determined that the features are coherent over long time periods so that we now show the results 
of taking the derivative on the non-smoothed instead of the smoothed fields. We have used the 
non-smoothed fields shown in figure ([6]) to take numerical derivatives of the Br trend lines. The 
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Fig. 7. — The time dependence of the radial component of the magnetic 
field Br for six selected latitudes near the sun's south pole. The dashed lines 
across the figure correspond to zero field for the top and bottom latitudes. 
Each latitude line is offset by 5 G from the adjacent line. 
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Fig. 8. — This figure gives a map of the time derivative of the radial magnetic field dB^/dt in units 
of G/yr. The color-bar on the right gives the coding of the mapped colors. The magnetic field for 
this mapping has only been averaged over a single Carrington rotation. Gaps in the record have 
been filled using spline interpolation with the filled portions being indicated as in figure Q. 
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Fig. 9. — This figure shows Br smoothed in the top panel and the difference 5Br between Br and 
the smoothed Br in the lower panel. The smoothing function is a truncated gaussian with a width 
of 2.5 years. Gaps filled by spline interpolation are indicated on both panels as in figure ^. 
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resulting time derivative map is shown in figure ([8]) . The rough regularity of the pattern in figure ^ 
was unexpected and represents a behavior that deserves attention. The continuity of the features 
from low to high latitudes is particularly striking and is found for even relatively weak ripples. 

The overall patterii of the ripples brought out by the time derivative is related to the structures 



found by IVecchio et al.l (|2012l ) who employed an iterative smoothing/differencing method called 
Empirical Mode Decomposition to isolate what are called Intrinsic Mode Functions (IMF's). The 
process includes spline fitting to maxima and minima of the starting function and a set of tests to 
determine when an Intrinsic Mode Function has been accurately found. The resulting functions are 
collected together in groups to show time and space dependence of the magnetic field. The choice of 
IMF's to group governs the range of frequencies retained. The authors show two groups, one with 
frequencies near one year and the other with frequencies between 1.5 and 4 years. Although both 
group results resemble our maps based on dBr/dt, there are differences that are significant such as 
the patterns that drift toward the equator on their maps - a drift which is not seen on our maps. 
We considered the use of spline fitting functions to remove the long-time background variations 
from our Bj. maps in order to display the ripples as magnetic field strength variations instead of 
using the time derivative but we found that the spline fit approach can introduce artifacts due to 
the stiffness of the fit function and prefer the method described in the following paragraph. 

We have adopted the difference 6Br between the raw Br (figure[6]) and a smoothed Bj. (figureO 
top panel) as a simple and effective way of representing the ripples. Figure [9] shows the result when 
we use a smoothing function with a truncated gaussian of width of 2.5 years. This choice of 
smoothing effectively divides the function into a primary wave, dipole-like component and a ripple 
component. At any fixed time, the SBj. function represents a perturbation on top of the primary 
component. The near-sun heliosphere is generally considered to be current f ree and describable 



by potential field solutions (JAltschuler fc Newkirklll969l : IWang &: Sheelevlll992l ). Consequently, the 



response of the field above the photosphere is linear to variations in the observed field. We can then 
consider the heliospheric field in this region to be the sum of three parts: the smoothed primary 
wave, a superposed response from the ripples and finally the longitude-dependent component that 
includes the effects of active regions and sunspots. The three vector components of the field include 
contributions from each of these components. If the ripples acted in isolation this response would 
take the form of loops whose foot points are separated by 15° to 30° of latitude and which would 
migrate mostly toward the poles. Compared to the smoothed wave, the amplitude of the loops is 
about 30% of the primary dipole-like field, at least during periods of high activity. However, the 
actual trace of the field lines is dominated by the more structured longitude-dependent component 
which is the strongest of the three. The effect of the ripples is to tilt the field lines so that they 
go in the direction of the loops that would be present if the longitude-dependent structures were 
missing. Detailed models are needed to determine if this effect is present. 

We searched for periodicities in the dBr/dt map using the Lomb-Scargle periodiogram ap- 



proach (IScargldll982l ). This method has the advantage that the time series does not need to be 



evenly sampled. Consequently, we are able to study the time series without gap filling and with- 
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Fig. 10. — This figure shows a sample of four periodograms selected to be for the latitudes 
between 40° and 60°N/S. The upper figures are for cycle 22 while the lower figures are 
for cycle 23. The identified frequencies are shown by the vertical dashed lines (red in the 
on-line version). 
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Fig. 11. — This figure shows the periods for all the peaks like those in figure (jlOp . The 
periods are sorted by the latitude band where they were found and differing plotting symbols 
are shown for each cycle. 
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out smoothing. The resulting periodograms for single latitudes are erratic with the peak location 
being function of latitude as well as a function of time when the original series is taken as several 
independent parts. A set of sharp peaks that are not stable is often an indication that the life- 
time of the underlying phenomenon is shorter than the full time series. In order to better match 
the time series to the lifetime of the magnetic pattern, we sub-divided the series into three parts 
corresponding roughly with the times of cycles 21, 22 and 23. In addition we have binned the 
latitudes into eight ranges, each of which is 20° wide. This gives 24 separate periodograms. There 
are few regularities among the periods found in these periodograms. A sample of four is shown 
in figure (|lUp . Each of the identified peaks has a vertical dashed line at the measured frequency. 
All of the periods for identified peaks are shown as symbols on figure ([TT]) grouped according to 
latitude band. The separate solar cycles are identified by the symbols according to the legend on 
the upper third of figure (jlip . There are no patterns among these periods and we conclude that we 
have not determined any organizing principle for the periodicity in the magnetic field fiuctuations. 



4. Discussion 

Maps of the sun's magnetic field like that of figure ([6]) are one of the most p owerful tools for the 



study of the solar cycle and have been used regularly since the publication by iHoward &: Labonte 



1981 ). Improved maps have beei i provided reg ularly with recent examples being found in lHathawa 



2010l ). IVecchio et al.l ( 20121 ) and iPetriel ( 20121 ). In most of these prior maps the patterns known as 



surges are clearly visible and usually lead directly to the reversal of the polar fields. Also visible 
in these is what can be called a counter-surge a few years after the main surge during which the 
polarity of the previous cycle has a return to dominance in the mid-latitudes (30°N/S to 60°N/S) 
but with inadequate strength to overcome the ne wly dominant po lar field. These counter-surges 



are most evident in our figure ([6]) and the map in iHathawayl ( 2010l ). 



Although on the full magnetic map the counter-surges appear to be relatively minor fiuctu- 
ations on the overall pattern of the solar cycle, they actually represent a fundamental problem. 
According to the standard picture of the cycle involves the sequence: a) the polar field generates a 
toroidal field from differential rotation, b) the toroidal field strengthens until it can generate bipolar 
sunspot groups, c) the bipolar sunspot groups are oriented according to Joy's Law such that their 
trailing spots are poleward of the leading spot so that the trailing spot polarity becomes dominant 
in the mid-latitudes by virtue of the spot's favorable position relative to the pole, d) the trailing 
spot polarity reverses the previous dipole field so that the sequence can repeat. As long as the Joy's 
Law tilt configuration of bipolar spots is stationary, the leading spot field cannot become dominant 
in the mid-latitudes independent of the number of spot groups present. During a counter-surge 
the leading spot field is dominant indicating that the spot tilt is opposite that given by Joy's Law. 
Although the counter-surges are rare and might be dismissed as anomalies, the ripples that are 
revealed by a close inspection of the magnetic field maps as enhanced in figures ([8]) and ([9]) show 
that in fact the polarity dominance in the mid-latitudes is continually varying. In the mid-latitudes 
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Fig. 12. — The map of the meridional circulation a s a function of time and latitude 
derived with the method described by lUlrichI (|201Cll ). Each 6-month bin is the result 
of averaging over the first or second half of each year. An even-odd feature present 
in the map was removed by using a (wj)mtcrod = 0.25t;j_i -|-0.50z;j-|-0.25fj+i temporal 
filter on each latitude. The final time slot corresponding to the last 6 months of 
2012 is in fact derived from only the months of July to Sept. 
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where there are few spots, our average over a Carrington rotation gives us a quantitative measure 
of the degree to which the traihng polarity is in fact doruinant . Although the tilt angle is, for a 
long-term average, a function of the latitude (JLi fc: Ulrichl 12012 ) the present study shows that the 
tilt angle is not steady. 

As is best seen in the time derivative plot of figure ([8]) the magnetic ripples alternate between 
reducing and enhancing the existing dipole structure with the cumulative effect causing the dipole 
component reversal. We can understand how the changes in the dominant magnetic field in mid- 
latitudes are relat ed to sunspot emergence by using the Babcock-Leighton model as extended by 
Wang et al.l (|l99lh . According to this model the time d erivative dBs/dt h as the form of a sum of 
an advection term, a diffusion term and a source term (jWang et al.lll99ll . see equations 5 and 6). 
The form of the source term, S, is: 

1 d /i?0a sin 7' 



S{RQ,9,t) oc 



(5) 



sin^a^ 

where Bg is the surface field strength which we identify with Br, B^j, is the toroidal field strength in 
the sub-surface zone where sunspots appear, a is the absolute value of pole separation between the 
components of a bipolar magnetic region, 7 is the Joy's Law tilt angle and r is the time scale for 
the features to erupt. The coefficient of proportionality is the ratio of several constant geometric 
factors. The dominant polarity in the active zone is generated by the source term from sunspot 
groups and this polarity then spreads poleward as a result of the advection and diffusion terms. A 
basic property of each sunspot cycle is the structure of the deep seated toroidal field B^ which does 
not change sign until the next cycle. Consequently, the alternating sign of dB^/dt can only be due 
to an alternating sign of 7. Thus the relative orientation of the spots as defined by the Joy's Law 
tilt angle 7 must be opposite from normal during periods of reversed sign of dBr/dt. 

The wave pattern of the magnetic ripples drifts in latitude as if their wave pattern is being 
advected by the meridional circulation. The actual ripple pattern is complex in many places and 
simple in a few. In those simple places the peaks and troughs can be tracked so that it is possible to 
determine a pattern drift velocity. While that pattern drift velo city is compa. rable in magnitude to 
the surface Doppler meridional circulation velocity reported bv lUlrichI (J20ld ). the match in detail 
i s qui te imperfect. For reference a 2-d map of the meridional circulation velocity found bv lUlrich 
([20101) is shown in figure [T2l The dual cell structure on this map represents a particular problem 
since it implies a zero drift rate along the boundary between the dominant poleward circulation cell 
and the polar cell where the circulation is equatorward. There are no cases where the drift of the 
magnetic ripple stagnates at this boundary. This disagreement indicates that the pattern drift is 
not caused by surface mass motions alone and is a consequence of a more complex interaction such 
as more deeply seated mass motion, magnetic forc es or flux emergence. An altern ate determination 
of the meridional circulation in a recent study by iRightmire-Upton et al.l (120121 ) based on pattern 
tracking of small-scale magnetic features found that the reverse cell at high latitudes was absent 
during the period from Apr. 2010 to Feb. 2011. 



A striking property of the ripples is their long duration - lifetimes of two to three years and 
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extending over latitudes starting at 20° and ending often at the pole. Some of the ripples start 
at the equator and extend into both hemispheres although this behavior is less common. The 
general time scale of the ripples is sir nilar to that of the Quasi-Bie i inial Oscillations that have 



been discussed by a variety o f authors (JHowe et al.ll200d : lBumball2003l : lKanell2005l : iJavaraiah et al 



20091 : iBroomhall et al.ll2012l ). The cause of these QBO's is not understood but the presence of a 
periodicity in this range in rotation rates, especially deep in the convection zone, suggests that the 
QBO's are related to deep-seated processes governing the solar cycle. 

It is interesting that in the northern hemisphere there are three wider than average bands of 
dBr/dt of -, + and - sign starting at about 60°N at times of 1978, 1988 and 1998 and ending near 
the poles at times of 1980, 1990 and 2000. These bands are of long enough temporal duration that 
they do not show up in the 6Br plot of figure [9] but rather are seen in the smooth B^ part of the 
figure and immediately preceed the reversal of the polar field. These correspond roughly to the 
times when coronal features unde rgo a rapid migration from mid-latitudes to the poles called the 
"rush to the poles" (jAltrocldl201ll ). The stronger-than-average ripples in the rate of magnetic field 



change are evidently associated with these coronal processes while the majority of the ripples do 
not have an impact on this global str ucture. I r i add ition, it is noteworthy that the rotation rate of 
the coronal features also discussed by lAltrockl (J2003l ) is significantly different from the rate derived 
from the Doppler shift in the polar regions. This is consistent with the idea that the magnetic field 
motions and evolution is not simply a consequence of advection by the matter on the solar surface. 



5. Conclusions 



We have presented improvements to the derivation of solar surface magnetic field strengths 
as measured at the Mt. Wilson Observatory. These improvements include the deduction of a tilt 
angle C for the orientation of the field lines near t he po l es and the filling of unseen portions of the 
poles using a method based on the techniques of iTranI (|2009l ). These improvements substantially 
reduce the annual variation in the line-of-sight field due to the changing viewing angle (the bo 
effect). Owing to the low-noise characteristics of the final magnetic field map, we are able to 
identify the presence of ubiquitous ripples in the magnetic field that generally start at low latitudes 
and propagate to the poles. The motion of these ripples cannot be due to advection at the surface 
alone due to the apparent inverse relationship between the surface meridional circulation and the 
apparent travel speed of the features. It is likely that in addition to the hydrodynamic influence of 
the convection zone matter, large-scale magnetic forces also play a role in the structure and motion 
of the ripples. 

The ripples in the dB^/dt map clearly are associated with the process that reverses the dipole 
during the solar cycle. Because we are considering averages over full circles at constant latitude, 
the Bj. value we find represents an axially symmetric pattern. At the time of the dipole reversal the 
value of dBy/dt is larger than average (note that the 5Br part of figure [9] does not show this effect 
because its duration is longer that 2.5 years and the deviation shows up on the Bj. smoothed panel 
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instead) and is associated with the "rush-to-the-poles" structural changes in the c orona. A part 



of dyn amo models of the solar cycle is the a-effect [see for example the review by ICharbonneau 



(J20ld )] whereby rising portions of the underlying toroidal magnetic field are twisted in such a way 
that the global dipole field is reversed. This twisting is manifested as the Joy's Law tilt of active 
regions and the magnetic ripple responds to changes in this Joy's Law tilt. 

It is possible that there are relati onships between th e magnetic ripples and other solar c ycle 



phenomena such as flare produ ction (|Rieger et al.lll984l ). sunspots (JKrivova fc Solankil |2002| ) or 



geophysical indices (JKand l2005l ). Since the ripple structure is large scale and long lasting, it 
could be associated with structures in the heliospheric magnetic field that are manifested in the 
geophysical indices. However, in the absence of a specific physical model as a basis for a search 
for such relationships, temporal concident and lagged cross-correlation studies may provide the 
best available tool. Because the magnetic ripples do not show any clear and stable periodicities, a 
coincidence of periods is only a rough guide to a search that would need confirmation through a 
cross-correlation approach. 
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